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Self-energy corrections due to fermion-magnon interaction are studied in the antiferromagnetic 
state of the t — t' — t" Hubbard model within the noncrossing approximation in the full U range from 
weak to strong coupling. The role of classical (mean-field) features of fermion and magnon dispersion, 
associated with finite U, t' , t" , are examined on quantum corrections to quasiparticle energy, weight, 
one-particle density of states etc. A finite-?/ induced classical dispersion term, absent in the t — J 
model, is found to play an important role in suppressing the quasiparticle weight for states near 
k — (0,0), as seen in cuprates. For intermediate U, the renormalized AF band gap is found to be 
nearly half of the classical value, and the weak coupling limit is quite non-trivial due to strongly 
suppressed magnon amplitude. For finite t' , the renormalized AF band gap is shown to vanish at 
a critical interaction strength Uc, yielding a spin fluctuation driven first-order AF insulator - PM 
metal transition. Quasiparticle dispersion evaluated with the same set of Hubbard model cuprate 
parameters, as obtained from a recent magnon spectrum fit, provides excellent agreement with 
ARPES data for Sr2Cu02Cl2. 



I. INTRODUCTION 

Driven by the intensive effort to understand the fasci- 
nating properties of doped cuprates, substantial theoret- 
ical progress has been made in recent years in the area of 
strongly correlated electron systems. In particular, the 
study of quantum spin fluctuations in the antiferromag- 
netic (AF) state and the renormalization of quasiparticle 
properties of doped carriers due to coupling with spin 
fluctuations has attracted much attention in the context 
of antiferromagnetic spin correlations, anomalous normal 
state properties, spin-fluctuation mediated pairing, opti- 
cal conductivity etc.^ 

Photoemission studies of quasiparticle properties of 
a doped hole in the antiferromagnet Sr2Cu02Cl2 have 
revealed^i^ nearly isotropic dispersion around k = 
(7r/2,7r/2), a bandwidth of around 0.3 eV, nearly same 
energy at (tt, 0) and (0, 0), and sharp quasiparticle peaks 
only near (7r/2,7r/2) and (27r/3, 0). In the context of 
photoemission studies, an overwhelming majority of the- 
oretical studies consider the t — J model, involving only 
nearest-neighbour spin coupling and excluding double oc- 
cupancy completely, thus retaining only the most signif- 
icant physics of the large U limit. Quasiparticle renor- 
malization has been obtained within the self-consistent 
Born (noncrossing) approximation involving scattering 
of mobile carriers due to multiple emission/ absorption 
of spin waves. Well defined quasiparticles are obtained 
near the top of the band at k = (7r/2,7r/2), with sig- 
nificantly reduced coherent weight, and a broad inco- 
herent background. The experimentally observed en- 
ergy difference between k = (7r/2,7r/2) and (tt, 0) states 
is too large to be understood within the t ~ J model. 
By including additional next-nearest-neighbour (NNN) 
and NNNN hopping terms t' and t", which control the 
(0, 0) - (tt, 0) and (0, 0), (tt, 0) - (7r/2, 7r/2) energy sepa- 
rations, respectively, substantially improved fittings have 
been obtained^'^iSiLSiSiiSiii with parameter values in the 



range: J = 0.12 - 0.15 eV, t = 0.3 - 0.4 eV, t' = 
-(0.07-0.14) eV, and t" = 0.05-0.09 eV. Quasiparticle 
dispersion fits have also been obtained within the t—t'—t" 
Hubbard model at the mean-field leveljiSii^ and within 
the Born approximation,^^ The need for more realistic 
microscopic models for cuprates which include additional 
hopping terms has also been acknowledged recently from 
band structure, photoemission, and neutron-scattering 
studies of high-Tc and related materials. ^^d^d'^d^d^i^" 
Estimates for \t' /t\ range from 0.15 to 0.5. 

While electronic excitations in cuprates can be quan- 
titatively understood within the t — J model with addi- 
tional hopping terms, certain features of the low-energy 
magnetic excitations (spin waves) cannot be understood, 
even qualitatively. Recent high-resolution inelastic neu- 
tron scattering studies of the spin-wave spectrum of the 
square-lattice spin-1/2 antiferromagnet La2Cu04 have 
revealed noticeable spin-wave dispersion along the AF 
zone boundary.^"'^ While the simplest explanation for the 
observed spin-wave dispersion is in terms of a ferromag- 
netic NNN spin coupling J', it was pointed ouliSiiS^ that 
a more natural explanation lies within the one-band Hub- 
bard model, for which the strong coupling expansion in 
powers of t/U generates extended-range spin couplings. 
Fits of the spin-wave spectrum yield t « 0.35 eV and 
U/t^7, with |f|/t = 0.25,aLSa 

Given that certain key features of spin-wave excita- 
tions in cuprates have now been ascribed to finite U, 
it would be reassuring if a unified description of both 
electronic and magnetic excitations in cuprates can be 
obtained with the same set of Hubbard model param- 
eters. The present work does provide such a unified 
description. In this paper we study the quasiparticle 
properties of a doped carrier in the AF state of the 
t — t' — t" Hubbard model within a mean-field-plus- 
fluctuation approach. Specifically, we study the self- 
consistent fermion self energy due to fermion-magnon in- 
teraction within the noncrossing (self-consistent Born) 
approximation. Not only is a Hubbard model approach 
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to quasiparticle renormalization evidently more realistic 
for cuprates, it also provides a complimentary perspec- 
tive involving a straightforward physical interpretation. 

The perturbative approach for the t — J model in- 
volves a slave-fermion representation to implement the 
single-occupancy constraint, mapping to a holon-spinon 
Hamiltonian, and then an expansion in the hopping 
term t. As the hole hopping is accompanied by up- 
turned or restored spins, the hole motion is renormalized 
due to its coupling with magnon emission and absorp- 
tion, and within the self-consistent Born approximation 
(SCBA), which is intimately related to the retraceable 
path approximation)^^ the hole self energy amounts to 
a self-consistent, second-order correction (in t) involving 
the hole-magnon scatteringi^^i^^i^^i^'^i^^i^^ Exact diago- 
nalization for small systems has also been carried outn^S 
On the other hand, the fluctuation approach adopted 
in this paper involves a many-body expansion in pow- 
ers of U, wherein the hole renormalization is physically 
due to the dynamical transverse fluctuations about the 
mean fleld, and the hole self energy amounts to a self- 
consistent, second-order correction (in U) involving in- 
elastic, spin-flip scattering induced by the transverse fluc- 
tuations. This fluctuation perspective of hole renor- 
malization is physically relevant for cuprate antiferro- 
magnets, where the observed Cu moment is reduced by 
nearly half due to strong spin fluctuations. Indeed, this 
fluctuation-induced quantum correction to the ordered 
moment is described by an interband self-energy con- 
tribution involving interband spectral weight transfer, 
which is absent in the t — J model approach. 

Reassuringly, the formally weak-coupling fluctuation 
approach for the Hubbard model, when carried to the 
strong coupling limit, yields (for the intraband self en- 
ergy) exactly the t — J model result of the strong cou- 
pling approach, thus continuously interpolating between 
the weak and strong coupling limits, and allowing for 
quasiparticle studies in the full U range. In the spin- 
density- wave (SDW) state for finite U, several important 
aspects of correlated electron systems, which are pro- 
jected out in the t — J model, can be addressed, such 
as interband spectral weight transfer, one-particle den- 
sity of states for both bands, renormalization of the AF 
band gap, and the possibility of vanishing band gap even 
at moderate U, thus providing an understanding of the 
antiferromagnetic insulator - paramagnetic metal tran- 
sition within a spin/charge fluctuation scenario. These 
finite-C/ aspects have not been systematically studied in 
the full U range in earlier studies. 

Earlier studies of self-energy and Green's function cor- 
rections due to spin fluctuations in the AF state of the 
Hubbard model have examined one-loop quantum correc- 
tions to sublattice magnetization, 3J*2SiSi2i2Si2& ground 
state energy;^ and quasiparticle dispersionj22i2i and 
later extended to the SCBA level with regard to quasi- 
particle dispersion and weight ^^s^ For intermediate and 
strong coupling, the one-loop results are only a poor ap- 
proximation, as they differ substantially from the self- 



consistent calculations.^^ Only the interband self energy 
contribution was considered in Ref. [38], whereas in- 
traband contribution was argued to be more important 
for t/Jy/S » 1 within a l/S* expansion in Ref. [14], 
where the SCBA self energy was studied in the con- 
text of photoemission results. A finite momentum cut- 
off in the magnon spectrum was introduced to model 
the suppressed weight of zone-boundary magnons, as 
indicated by two-magnon Raman scattering. As the 
high-energy, zone-boundary magnons contribute signif- 
icantly to the SCBA self energy, their elimination biases 
the quasiparticle dispersion towards the mean-field re- 
sult, which has been noted to yield, for t' = — 0.5J, 
a surprisingly good description of the photoemission 
results Good agreement with experimental data 
was obtained for J/t = 0.4 and t' — —0.4 J, with a 
quasiparticle bandwidth of around 3 J. Longitudinal and 
charge fluctuations are negligible in the strong coupling 
limit pii^^i^ and all fluctuations are small in the weak 
coupling limit. 

Fluctuation self energy has also been examined in the 
normal metallic state, close to the transition to the AF 
insulator, by assuming a model dynamical spin suscep- 
tibility x(g, cj),'^- In this case, instead of the intra and 
interband self energies, one has intra and inter Fermi 
surface contributions from scattering states k — q lying 
on same or opposite side of the Fermi surface. With 
increasingly long-range AF correlations, the dynamical 
spin susceptibility is more sharply peaked at q = Q, the 
AF wavevector, and the inter Fermi surface contribution 
dominates as wavevectors k and k — Q lie on opposite 
sides of the Fermi surface. Formally resembling second- 
order correction, the inter Fermi surface self energy shifts 
the particle (hole) energy upward (downward), leading to 
a pseudo gap formation in the electronic spectrum with 
increasing U. The pseudo gap develops into the AF gap 
as long range spin order sets in, thus providing a compli- 
mentary approach from the paramagnetic side. 

Both SCBA studiesi^iS assume the strong coupling 
structure of the magnon dispersion and amplitudes, 
whereas in the weak coupling limit, magnon modes with 
energy approaching the band gap have strongly sup- 
pressed amplitudes. The flnite-J7 aspects mentioned 
above, such as DOS and AF band gap, have therefore 
not been systematically studied in the full U range. 

In this paper we examine the finite-?/ aspects in terms 
of an interplay between classical (mean-field) features 
of fermion and magnon dispersion and quantum (self- 
energy) corrections on quasiparticle dispersion. The in- 
terplay between classical and quantum effects is man- 
ifested in several ways. For example, the AF state 
of the Hubbard model at strong coupling involves a 
mean- field dispersion term 4J7^, where J — AP/U and 
7k — (cos kx + cos ky)/2. Absent in the t — J model, this 
large classical dispersion term has an energy scale twice 
that of the magnon energy, and therefore quasiparticle 
renormalizations in the Hubbard and t — J models are 
expected to differ with increasing J. Furthermore, the 
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finite- ?7 induced dispersion term 4J7^ intrinsically con- 
tains effective NNN {f = J/2) and NNNN {t" = J/4) 
hopping terms, and these finite-C/ contributions are en- 
tangled with the physical hopping strengths in the t' and 
t" values determined by quasiparticle dispersion fitting 
for cuprates within the t — t' — t" — J model. 

Another manifestation of the interplay is a competition 
between classical and quantum effects. Thus, while quan- 
tum correction lowers the hole energy near (7r/2,7r/2) 
in the t — J and Hubbard models, the dispersion term 
—At' cos kj. cos ky favours the (tt, 0) state (for positive t'). 
An effective one dimensional dispersion can thus result 
from a net cancellation of the dispersion along the mag- 
netic zone boundary. 

An instance where both classical and quantum effects 
come together involves the AF band gap for particle- 
hole excitations. Both, the classical dispersion term 
—At' cos kx cos ky and the quantum correction, reduce the 
band gap significantly, which may therefore be expected 
to vanish at a moderate U value. The antiferromagnetic 
insulator - paramagnetic metal transition associated with 
the vanishing band gap, relevant for transition-metal ox- 
ides such as V2O3, can thus be explored within a spin- 
fluctuation theory. 

Competition between the two classical dispersion 
terms in the t — t' Hubbard model introduces asymme- 
try in the two AF bands, in terms of bandwidth and 
density of states. The self-energy corrections for added 
holes and electrons are therefore expected to be signifi- 
cantly different. Indeed, we find a dramatically enhanced 
strong-coupling signature in the self-energy correction for 
the narrow band. 

After a brief review of the mean-field AF state of the 
t — t' — t" Hubbard model (section II) , the intraband self 
energy for one added hole (electron) is obtained in the 
noncrossing approximation (section III). Various quasi- 
particle properties are evaluated for the Hubbard model 
in the full U range (section IV). The role of finite t' and 
t" is studied in section V, and our conclusions are pre- 
sented in section VI. Some general features associated 
with the self energy in the AF state are put in Appendix 
A, and the interband self energy contribution is discussed 
in Appendix B. 



II. MEAN-FIELD AF STATE OF THE t - t' - t" 
HUBBARD MODEL 

To introduce the notation, we briefly review the mean- 
field (Hartree-Fock) state of the Hubbard model 



Hartree-Fock (HF) Hamiltonian 
H^{k) = (ej, + e'^)l + 



(2) 



for spin cr, in terms of the free-particle energies ek = 
—2t{coskx + cos ky), ej^ — — At' cos kx cos ky, and e'^ — 
—2t"{cos2kx + cos2ky). Here 2A — mil, where m is the 
sublattice magnetization. 

For the classical level (mean-field) fermion propagator. 
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for spin t, and [Gj'(k,a;)] — [iTa;jL^| 
sublattice symmetry, where [a^] = 
AF band energies 
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While the 



(4) 



for the lower and upper bands are modified by t' and t" , 
the fermion amplitudes |k) = (ctk Pu.) remain unchanged: 



al 




(5) 



In the strong coupling limit [U >> t), the majority and 
minority densities reduce to 1 — 4/^^ ^^'^ ^k/^^- 
Similarly, for the magnon propagator we have 



UqVq 



q 

u„v, 



UnV. 



q^q 



n-rir, 



IT] 



O + flq — irj 



(6) 



at the classical (ladder-sum) level. For the Hubbard 
model in the strong coupling limit the magnon ampli- 
tudes |q) = (wq Vq) and energy fiq are given hy^ 



q 



(Vy^ + i) /2 



l/v/l-7^-l)/2 



(-7q/yr^) /2 

2J, 



where 7q 



'1 
(cos qx 



7q 

- COS qy)/2 



(7) 



(1) 



on a square lattice, with hopping terms ts — t, t' , and t" , 
connecting NN, NNN, and NNNN pairs of sites {i,i + 6), 
respectively. It is convenient to use the two-sublattice 
basis as translational symmetry is preserved, and the 



HI. SELF ENERGY IN THE NONCROSSING 
APPROXIMATION 

We now obtain the self energy due to fermion-magnon 
interaction in the noncrossing (self-consistent Born) ap- 
proximation (Fig. 1), in which all noncrossing diagrams 



4 



I(k,CO) 




FIG. 1: Self energy in the noncrossing approximation. Wavy 
lines represent the magnon propagator. The bare fermion- 
magnon interaction vertex is U. 



are summed, representing multiple magnon emission and 
absorption processes. The magnon lines are consid- 
ered here at the classical (ladder-sum) level. The first- 
order quantum corrections to magnon propagator yield 
momentum-independent renormalizations of the magnon 
amplitude and energy in the strong coupling limit, "^'^ 
which can be incorporated easily. 
The SCBA self-energy 

[E(k,c.)] = [/2^ / |^[x-+(q,f^)][Gi(k-q,c.-l})] 

(8) 

for state |k |) involves the renormalized propagator for 
state |k — q j). However, from spin-sublattice symmetry, 
the fcrmion propagators are identical if their spin and 
sublattice indicies are simultaneously switched, so that 



[Gi(k-q,c.-0)] = 



^ [gx]|k-q) (k-q\[a^] 

iu-n- £;o_q - Sk-q(cj - n) ' 

(9) 

using Eq. (Al). The above replacement in Eq. (8) thus 
yields a self-consistent equation for the fermion self en- 
ergy Sk('j-') which is independent of spin. In the fol- 
lowing, we consider only the dominant intraband self en- 
ergy contribution in which the intermediate fermion state 
|k — q 1) lies in the same band as the state |k t). The 
intcrband contribution is discussed in Appendix B. Due 
to an exact cancellation with the Hartree self energy, the 
intcrband self energy does not contribute to gap renor- 
malization in the strong coupling limit, and the quasi- 
particle coherent weight reduction is small, even in two 
dimensions. 



A. One hole (electron) in the AF 

The situation is simplest for a single hole introduced in 
a filled lower AF band, as the available scattering states 



|k — q) span all of the lower AF band. From particle- hole 
symmetry, the situation is symmetrical for an added elec- 
tron to the empty upper band, under sign change of t' 
and t" . For finite doping, there are two intraband self en- 
ergy contributions arising from scattering states |k — q) 
above or below the Fermi energy, involving retarded or 
advanced magnon poles in Eq. (8) The local reduction in 
magnetic order near the hole is due to intraband spectral 
weight transfer across the Fermi energy. 

The motion of a hole in an antiferromagnetic back- 
ground has been studied extensively in the recent pasti 
The string of upturned spins and broken antiferromag- 
netic bonds left in its wake renders the hole motion highly 
incoherent, and the description of the hole carrying along 
the locally scrambled AF spin arrangement is a highly 
nontrivial many-body theoretical problem. 

Several aspects of the hole propagator, such as 
the spectral function, quasiparticle dispersion, band- 
width and weight, have been studied for the t — J 
model within the noncrossing (self-consistent Born) 
approximationi^^iS*2SiSL22iSS Comparison with exact di- 
agonalization resulta^S show the SCBA results to be in 
good agreement for small J, and contributions of crossing 
diagrams have been argued to be smallii^iSi 

As Gx(k — q, w — n) in Eq. (8) is an advanced propa- 
gator (for the intraband contribution) , representing scat- 
tering states for the added hole, the Vl integral (replaced 
by a contour integral over the upper half plane) picks a 
contribution only from the advanced part of the magnon 
propagator in Eq. (6). It is convenient to introduce a 
composite amplitude 

|a(k,q)) = |q)x[a.]|k-q)=(^;;^fj^-;j) , (10) 

for the magnon-hole interaction vertex, in terms of which 
we obtain for the hole self energy 



Ek(c^) ^ (k|[S(k,c.)]|k) = (a^ /3£)[S(k,c.)] 



/3k 



U 



q 



(k|a(k,q))(a(k,q)|k) 



(«kMq /3k-q + A^'q ^k-q)^ 
f^q - -^k-q - Sk-q(cc' + r!q) 



(11) 

The above yields a self-consistent equation for the hole 
self energy in the form of an integral equation, which de- 
scribes the hole renormalization due to multiple magnon 
emission and absorption processes as it moves in the AF 
background, yielding a quantum correction to the clas- 
sical (mean-field) hole dispersion. While this intraband 
self energy significantly redistributes spectral weight in 
the hole (particle) band, leading to bandgap renormaliza- 
tion, it leaves the integrated spectral weight unchanged, 
so that the spin densities and ordered moment are unaf- 
fected. Within the noncrossing approximation, Eq. (11) 
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is valid for arbitrary [/, in the whole range from weak to 
strong coupling, and for extended-range hopping terms. 
As the self energy depends on the classical dispersion, ex- 
plicitly through ^'k-q ^'^^ implicitly through the magnon 
energy Qq, an interesting interplay is expected between 
classical and quantum effects on quasiparticle energies. 



IV. SELF ENERGY FOR THE HUBBARD 
MODEL 

A. Strong coupling limit 

We now consider the self energy for the Hubbard model 
in the analytically simple strong coupling limit (U/t — > 
oo), and compare with the t — J model result. We write 
the classical (mean-field) hole energy 



A2 



-A- 



4Jc7k 



(12) 



where a classical energy scale Jc — 4:t'^/U has been intro- 
duced to distinguish the classical dispersion term 4Jc7^ 
generated at strong coupling. Substituting in Eq. (11) 
the fermion and magnon amplitudes and energies from 
Eqs. (5), (7), and (12), and shifting the energy by A to 
bring the zero of the energy scale to the lower band edge, 
the self energy is compactly written as 



(Mq7k-q + Vqlkf 



4>/c7Lq - Sk-q(^-' + f^q) ' 



Sk(c.)=tVV- 
— ' L 

q 

(13) 

which is exactly the SCBA result for the t — J 
modelr^'2£'2& apart from a classical dispersion term ab- 
sent in the t — J model. This shows that the formally 
weak-coupling fluctuation approach, involving a many- 
body expansion in powers of U, when carried to the 
strong coupling limit, yields the strong-coupling result for 
the hole self energy as well. Similar exact correspondence 
has been noted earlier for magnon dispersion, both at the 
classical level and including quantum corrections, sub- 
lattice magnetization, ground-state energy, two-magnon 
Raman scattering etc. 

While the effective fermion-magnon scattering matrix 
element is reduced from 0{U) to 0{t) due to the minority 
fermion amplitude /3, the gapless energy denominator in 
this intraband term leads to a self-consistent self energy 
of 0(t), which is much greater than the residual inter- 
band contribution of 0(J), surviving the leading-order 
cancellation (see Appendix B). 



B. Strong-coupling results 

The self-consistent numerical evaluation of the self 
energy was carried out on a 52 x 52 grid in k space, 
and a frequency interval Alo = 0.1 for w in the range 
— 10 < uj < 10. The self energy was iteratively evaluated, 
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FIG. 2: Variation of self energy with iterations, effectively 
illustrating the role of successively higher-order magnon pro- 
cesses. Here J = 0.5 and k — (7r/2, 7r/2). 




FIG. 3: The hole spectral function ylk(ti') shows the quasi- 
particle peak at the intersection of tj — with Re T,k{uj). 



starting with Sk(^) = 0. Typically, self-consistency was 
achieved within ten iterations. For Jc = 0, our numeri- 
cal results are in exact agreement with the t — J model 
results of Liu and Manousakisi2i 

The hole self-energy variation with successive itera- 
tions (Fig. 2) effectively illustrates the role of multiple 
magnon emission and absorption processes. The first it- 
eration yields the single magnon contribution, and suc- 
cessive iterations incorporate higher-order magnon pro- 
cesses. The significant difference, especially for J « t, 
between the first-order and the self-consistent results, 
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FIG. 4: Quasiparticle dispersion along different symmetry di- 
rections, for J = 0.5. 
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FIG. 5: Quasiparticle dispersion along the F — X direction 
with (+) and without (x) the classical dispersion term Jc, 
along with the mean-field dispersion. 



indicates the relative importance of multi-magnon pro- 
cesses. 

Prom Eq. (Al), intersection of w — with Re 
yields the renornialized quasiparticle energy where 
the spectral function peaks, as shown in Fig. 3 for k = 
(^,^). The quadratic growth of Im Ek(a;) near u) « 
E^ is due to the long-wavelength magnon modes. In 
Eq. (13), neglecting the weak dependence of 7k;-q and 
Sk-q in comparison with the linear dependence Oq = 
V^Jq of the magnon energy for small q, one obtains 

ImEkM ~ .V/,ci,(l/,)fe+,,m. + l.q + ReEk) 

- {^ + Etf. (14) 

Magnon damping at finite doping due to intraband 
particle-hole excitations is therefore expected to be im- 
portant for quasiparticle damping. 

The large classical term ^Jdt suppresses the self- 
energy correction in Eq. (13), in addition to provid- 
ing dispersion at the mean-field level, and these com- 
bined effects of Jc on quasiparticle dispersion are dis- 
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FIG. 6: The k = (0, 0) spectral function is significantly mod- 
ified by the finite-[/ induced classical dispersion term Jc, as 
shown in terms of the intersection of Re Ek(a)) with the line 
-I- 4Jc7k. Here J = 0.5. 



cussed below. Figure 4 shows the quasiparticle disper- 
sion along different symmetry directions for J = 0.5. 
Here the quasiparticle energies are obtained from the 
peak in the spectral function ^k(w) occur ing at lowest 
hole energy. The small dispersion developed along the 
M — M' direction (magnetic zone boundary) is purely 
quantum mechanical. The substantially different curva- 
tures at k = (7r/2, 7r/2) along the T - X and M - M' 
directions shows a highly anisotropic quasiparticle dis- 
persion with different effective masses. The quasiparticle 
energy is maximum for k = (7r/2,7r/2). 

Figure 5 shows a comparison of the quasiparticle dis- 
persion along the F — X direction with and without Jc, 
providing a comparison with the SCBA results for the 
t — J model (Jc = 0). The classical dispersion term 
Jc reduces the quasiparticle energy renormalization, and 
significantly flattens the dispersion near the F and X 
points, and also increases the quasiparticle bandwidth. 
Furthermore, for k = (0,0), the spectral weight of the 
highest-energy peak (lowest-energy for hole) is drastically 
reduced by Jc, in sharp contrast with the Jc = case, 
as shown in Fig. 6. This finite-[/ induced feature is 
quite significant with regard to photoemission studies of 
cuprates, where the k = (0,0) peak is seen to be ex- 
tremely weak in comparison with (7r/2,7r/2). 

The enhancement in the self-energy correction with de- 
creasing J is seen in Fig. 7, both through (a) decreas- 
ing quasiparticle weight for k = (7r/2, 7r/2) signifying en- 
hancement in incoherent part due to multi-magnon pro- 
cesses, and (b) increasing quasiparticle energy renormal- 
ization. Comparison with the t— J model results {Jc = 0) 
illustrates the increasing suppression of quantum correc- 
tion by the classical dispersion term with J. 
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FIG. 7: Comparison of quasiparticle weight and energy 
E^. with and without Jc, showing increasing differences with 
J. 
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FIG. 8; Comparison of strong-coupling and finite-17 spectral 
functions for k = (7r/2,7r/2) and (0, 0). 



C. Intermediate and weak coupling 

Self energy for finite U is of special interest as sev- 
eral additional features of the AF state can be explored. 
These include one-particle density of states for both 
bands, renormalized band gap, and the possibility of van- 
ishing band gap and a metal-insulator transition driven 
by spin fiuctuations. Furthermore, the magnon nature 
is qualitatively modified at finite f/, as studied earlier in 
detail in the intermediate)^ and weak^ coupling limits. 
Extended-range spin couplings generated at intermediate 
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FIG. 9: Increasingly incoherent hole spectral function with 
increasing U, for k = (7r/2,7r/2). 




FIG. 10: Approximate scaling of the real part of self en- 
ergy, showing a dominant single-magnon contribution with 
decreasing U. 



U remove the degeneracy in the magnon spectrum along 
the AF zone boundary. Also, the decreasing AF band 
gap 2A begins to play a significant role, and for small U 
the magnon amplitude is strongly suppressed for modes 
with energy approaching 2A.'*° These features were not 
explored in earlier SDW-SCBA calculationsi^iS for the 
Hubbard model, where only strong coupling expressions 
for magnon amplitude and energy were employed, and are 
naturally beyond the scope of the t — J model in which 
the upper AF band is projected out, and the magnon 
dispersion corresponds to NN spin coupling only. 

For arbitrary C/, the hole self energy SkC"^) is deter- 
mined from Eq. (11), with classical fermion amplitudes 
and energies taken from Eqs. (4) and (5), and the classi- 
cal (ladder-sum) magnon energy fJq and amplitudes Mq, 
Wq determined numerically, as described earlier 

Comparison of spectral functions evaluated for inter- 
mediate U with strong coupling results (Fig. 8) shows 
a noticeable reduction in the quasiparticle bandwidth. 
With decreasing [/, the hole spectral function becomes 
increasingly coherent (Fig. 9), the self energy approxi- 



8 




0.4 - [7 = 2 




-10 -8 -6 -4 -2 2 4 6 8 10 



1 1 1 1 

U = 4 


1 


1 1 1 1 


1 -j/^^'^'^ \ 


V 





-10 -8 -6 -4 -2 2 4 6 8 10 



0.5 I 1 1 1 ni 1 in 1 r 



0.4 - V = % 




-10 -8 -6 -4 -2 2 4 6 8 10 



FIG. 12: Renormalized one-particle (tunneling) density of 
states for U = 8, 4, 2. The corresponding mean-field gaps 
are 2A = 7.2,2.8,0.75. 



mately scales as (Fig. 10), and the shape of Rel]k(^^) 
increasingly resembles the single-magnon result (Fig. 3), 
indicating dominant lowest-order (J7^) contribution, as 
expected. Evolution of the k = (0, 0) spectral function 
(Fig. 11) shows that the peak at lowest hole energy re- 
mains nearly incoherent down to very low \J values. 

The one-particle (tunneling) density of states is shown 
in Fig. 12 for IJ = 8,4,2 (corresponding to mean-field 
gaps 2 A = 7.2, 2.8, 0.75), displaying the smooth approach 
towards the free-particle DOS with decreasing U . The 
band gap remains finite no matter how small U is. States 
are pulled within the classical energy gap 2A due to 
the fermion-magnon scattering, and the resulting fermion 
spectral function incoherence spreads the DOS over a 
broader frequency range. Pseudo-gap feature appears in 
the strong coupling limit. These DOS features are in 
good agreement with exact diagonalization results. "^^ 

For all [/, the state k = (7r/2,7r/2) continues to be 
the lowest-energy hole state, and therefore, from particle- 
hole symmetry, also the lowest-energy electron state in 
the upper band. The minimum electron-hole pair exci- 
tation energy therefore yields the renormalized band gap 
= 2|£;;j|, with k = (7r/2,7r/2). We find that the ratio 
E*/E^ is reduced to nearly 1/2 in a broad U range, as 
shown in Fig. 13, and approaches 1 in both weak and 
strong coupling limits, although with fundamentally dif- 
ferent behaviour. The weak-coupling trend indicates that 
the band gap remains finite for all U. For U > 8, we have 
used the strong-coupling self energy from Eq. (13). 

V. ROLE OF HOPPING TERMS " 

Finite hopping terms t' , t" add an additional dimension 
to the competition between classical and quantum effects 
on quasiparticle dispersion. Thus, the classical dispersion 
term e^. introduces a dispersion along the MM' direction 
of the Brillouin zone, renders the two AF bands asym- 
metrical, yielding narrow and wide bands for a given U 
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FIG. 14: Classical (a) and renormalized (b) one-particle den- 
sity of states for one added hole and electron, showing the 
vanishing of the energy gap. The classical- level band asym- 
metry strongly influences the quantum corrections. 



value, and also significantly reduces the band gap, mak- 
ing it possible for quantum fluctuation effects to close 
the gap at a moderate U value. Furthermore, the frus- 
tration due to competing AF spin couplings generated 
by t',t" leads to magnon softening, and the reduced 
magnon energy scale is expected to enhance the self- 
energy correction. This frustration effect is absent within 
the t — t' — t"—J model having only NN AF spin coupling. 

We have evaluated the hole self energy for finite U 
and positive t' using Eq. (11), with appropriate classi- 
cal fermion energies and amplitudes given in Eqs. (4) 
and (5). The quasiparticle energy difference between the 
{tt/2, 7r/2) and (tt, 0) states decreases nearly linearly with 
t' , and vanishes for t' « 0.1, resulting in a nearly one 
dimensional quasiparticle dispersion along the F — X di- 
rection. For higher t' values, the classical effect dom- 
inates, and k = (tt, 0) is the lowest-energy hole state. 
From particle-hole symmetry, these results also apply to 
an added electron in the upper band for negative t' . 

Figure 14 shows the one-particle density of states for 
t' — 0.5 and J7 « 5. While a substantial classical gap re- 
mains, the renormalized band gap is seen to just vanish. 
However, as the lowest-energy hole and electron states 
correspond to different momenta — (tt, 0) for hole and 
(7r/2, n/2) for electron — the band gap is indirect, yield- 
ing a finite optical gap. A further reduction in U will 
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FIG. 15: Spectral function shows sharp quasiparticle peaks 
for k — (7r/2,7r/2) and (27r/3, 0), as seen in photoemission 
experiments on cuprates. 
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FIG. 16: Quasiparticle dispersion along (tt, 0) — (0, tt) with 
same set of parameters as in Fig. 15, along with experimental 
ARPES data for Sr2Cu02Cl2 from Ref. [2]. 



lead to gapless particle-hole excitations with same mo- 
mentum. 

For one added electron to the narrow upper band, the 
characteristic strong-coupling signature of multi-magnon 
processes and the large energy shift relative to the mean- 
field DOS confirms the strong fermion-magnon scattering 
associated with a narrow band. The pseudo-gap feature 
is also significantly enhanced. On the other hand, for one 
added hole to the wide lower band, significantly weaker 
fermion-magnon scattering is evident. 

Taking the same cuprate parameters as obtained from 
a recent magnon spectrum fit,^^ and including a small 
t" term, we have obtained the spectral function for key 
states and the quasiparticle dispersion along (tt, 0) — 
(0,7r), as shown in Figs. 15 and 16. Excellent agreement 
with the photoemission experiments is obtained with re- 
gard to i) quasiparticle dispersion along (7r,0) — (0,7r), 
ii) sharp peaks near (7r/2,7r/2) and (27r/3,0), and iii) 
insignificant (low) spectral weights for (0, 0) and (tt, 0) 
states. With the same NN hopping term t = 0.35 eV, 
our finite-f/ study yields f « -0.09 eV and t" « 0.05 
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eV, which are quite smaller than the t — J model values 
(f « -0.12 eV and t" « 0.09 eV) of Ref. [3]. 



VI. CONCLUSIONS 

We have obtained one-particle self-energy corrections 
in the AF state of the t — t' — t" Hubbard model within 
a physically transparent fluctuation approach involving 
the dynamical transverse spin fluctuations. The full self 
energy consisting of the inter and intraband contribu- 
tions, describes both the quantum correction to ordered 
moment and ground state energy due to interband spec- 
tral weight transfer, as well as the intraband, inelastic 
spin-flip scattering leading to quasiparticle renormaliza- 
tion due to multimagnon emission and absorption pro- 
cesses associated with hopping of the mobile carrier. 

In the strong coupling limit, the intraband self-energy 
expression is shown to reduce exactly to the SCBA re- 
sult for the t — J model, apart from a classical dispersion 
term absent in the t — J model. That the formally weak- 
coupling many-body scheme can be carried over to the 
strong-coupling limit, indicates its applicability in the full 
U range from weak to strong coupling. By comparing the 
Hubbard and t — J model results, we have shown that 
the classical dispersion term has a significant effect on 
quasiparticle energy, mass, and weight, especially with 
increasing J, and indeed, some features of photoemis- 
sion experiments on cuprates, such as vanishing spectral 
weight for k = (0, 0), find natural explanation within the 
Hubbard model results. 

The Hubbard model study provides for a smooth in- 
terpolation between the strong and weak coupling limits. 
The importance of multi-magnon processes and string 
excitations in the antiferromagnet, manifested in pro- 
nounced oscillations in spectral function and density of 
states, gradually diminishes with decreasing [/, and the 
spectral function becomes increasingly coherent, the real 
part of self energy shows a dominant single- magnon con- 
tribution, and the DOS smoothly approaches the free- 
particle limit. For the Hubbard model, the one-particle 
density of states shows a finite band gap for all U values, 
and the ratio E* / E'^ of the renormalized band gap and 
the classical (mean-field) gap is found to be nearly 1/2 for 
a broad range of U values in the intermediate coupling 
regime. 

The role of fermion scattering states on quasiparti- 
cle renormalization due to fermion-magnon interaction 
is highlighted by the f'-induced asymmetry, leading to 
quite different self-energy corrections for the two AF 
bands. We find pronounced strong-coupling signature 
of multi-magnon processes for the narrow band, whereas 
the broad band exhibits significantly weaker renormal- 
ization. 

For finite i', we find that the AF band gap vanishes at 
a critical interaction strength Uc- For U < Uc, electron- 
hole pairs can be spontaneously excited due to band over- 
lap, which will reduce the sublattice magnetization and 



hence the classical gap, and thus further increase the 
band overlap. A first-order AF insulator - PM metal 
transition is therefore expected a,tU = Uc- The depen- 
dence of self energy suggests a significantly stronger sup- 
pression of the AF band gap in three dimensions, yielding 
a moderate Uc for even smaller t' value. 

With the same set of cuprate parameters as obtained 
from a recent magnon spectrum fitfSS excellent agreement 
with ARPES data for Sr2Cu02Cl2 is obtained, both with 
respect to quasiparticle weight as well as dispersion, ef- 
fectively providing a unified description of magnetic and 
electronic excitations in cuprates. The i" value required 
is only half of that obtained within the t — J model study 
of Ref. [3], highlighting the role of the finite-/7 induced 
classical dispersion term ^J"i^- 

In a three dimensional system at finite temperature, 
thermal excitation of magnons will enhance the self- 
energy correction, further decreasing the band gap and 
leading to, for U ^ Uc, a temperature driven first-order 
phase transition. An exploration of both quantum and 
thermal phase transitions from the AF insulator to the 
PM metal, as relevant for systems such as V2O3, is thus 
possible from quasiparticle renormalization due to self- 
energy correction. The t—t' Hubbard model also provides 
a concrete realization of a stable antiferromagnetic state 
for finite doping. The doped AF state is character- 
ized by magnon damping due to decay into particle-hole 
excitations across the Fermi energy^ which will have 
a significant role on quasiparticle damping through the 
imaginary part of the fermion self energy. These conse- 
quences of dimensionality, finite temperature, and doping 
are presently under study. 



APPENDIX A: SELF ENERGY IN THE AF 
STATE 

If [S(k,w)] represents the self-energy matrix in the 
two-sublattice basis, then in terms of the mean-field 



eigenvalue E^ and eigenvector |k) 
malized Green's function 

[G(k,c^)] 



/3k 



the renor- 
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[GO(k,( 



[S(k,c.)] 
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where expectation value of the self-energy matrix 



I]kM^(k|S(k,c.)|k) = (a;; PI) 



Saa Sab 
Sba ^bb 



(Al) 



/3k 
(A2) 



yields the self energy Ek(w) for state k. 

Quasiparticle dispersion, weight, and density of states 
can then be obtained from Eq. (Al) for the renormal- 
ized fermion propagator. The renormalized quasiparticle 
energy E^ is obtained by solving 
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Re Sk(w) , 
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FIG. 17: First-order self-energy corrections, representing (a) 
correction to the mean-field (Hartree) energy and (b) mixing 
between states of the two AF bands. 



and the quasiparticle weight is given by 

-1 



1 - ^Re I]k(w) 
ouj 



(A4) 



is the expectation value of the total self-energy matrix in 

the upper band state |k — q) = ^ ^- Therefore, 

for the lower band (hole) state |k) in consideration, the 
additional interband self energy 

EUc.) ^ (k|[E'(k,^)]|k) 



U 



q 



^ - ^^q - - K-^{U^ - f^q) 



(B3) 



The above interband contribution transfers spectral 
weight from lower to upper band, and together with Eq. 
(11), constitutes the full hole self energy, which describes 
both the hole renormalization due to multi-magnon in- 
traband scattering described earlier, as well as the inter- 
band transfer of spectral weight responsible for density 
and magnetization quantum corrections. 



In the absence of an intersection, the quasiparticle energy 
is obtained from the location of the peak in the spectral 
function 

A^iiu) = -Im Tr[G(k,c^)] = ^Im f ^ ' 

(A5) 

which also yields the density of states 

N{lu) = i y Im Tr[G(k, uj)] = V Ak(t^) . (A6) 

k k 

APPENDIX B: INTERBAND CONTRIBUTION 

An additional interband contribution to the self energy 
is obtained when the intermediate fermion state k — q in 
Eq. (8) lies in the opposite band. Absent within the 
t — J model approach, this contribution involves inter- 
band transfer of spectral weight, resulting in sublattice 
magnetization reduction due to quantum spin fluctua- 
tions. At the one-loop level, this interband contribution 
is exactly cancelled to leading order by the Hartree self 
energy correction due to spin density (ordered moment) 
renormalization. We also obtain the one-loop correction 
to ground state energy as a simple application. 

With upper band fermion state |k — q I), and a re- 
tarded magnon pole from Eq. (6), the interband contri- 
bution to the self-energy matrix of Eq. (8) is 



[E'(k,c.)] = C/^^ 



|a'(k,q))(a'(k,q)| 



(Bl) 

where, in analogy with Eq. (10), the composite fermion- 

/ u (Xif^ \ 

magnon amplitude |o-'(k, q)) — [ o "'^ ) , and 

Y Itq Pk-q J 

- f^q) = (k - q| [E(k - q, - |k - q) (B2) 



1. One-loop correction — Strong coupling limit 

The one-loop self-energy diagrams shown in Fig. 17 
represent (a) static correction to the mean-field potential 
due to density corrections, and (b) mixing between states 
of the two AF bands. Green's function correction due 
to diagram (b) yields one-loop quantum corrections to 
spin densities and sublattice magnetizationi^SiSi Spin-| 
density increases (decreases) on A (B) sublattice sites due 
to transfer of spectral weight from (to) the upper band. 
For diagram (a), these density corrections of magnitude 
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1 I yield a Hartree self-energy 
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(B4) 

And for diagram (b) representing mixing, for a state 
|k 1) in the lower (upper) band, and the intermediate 
state |k— q |) in the upper (lower) band, substituting 
the strong coupling fermion and magnon ampitudes in 
Eq. (Bl), we obtain the self energy for the lower band 



[Ef)(k,c.)P = C/2l^ 



1 



2^ UJ-Qq- E^^^ + IT] 



1 _ k-q 
J- U2 



- 1 
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and for the upper band 



k-q 



12 



-7q 
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"k-q 

(72 
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(B6) 



These self-energy corrections have the formal structure of 
a second-order process involving mixing between states, 
and result in energy lowering (raising) for lower (upper) 
band states. 

An exact cancellation to leading order (U) is seen be- 
tween the Hartree shift [sj"''] of Eq. (B5) and the mixing 
contributions [sj^''] of Eqs. (B6) and (B7), when the self 

energies T,^{uj) and S^(w) are evaluated from the corre- 
sponding matrices for the lower and upper bands, with 
uj — and , respectively, and the surviving terms 
of order J yield (for the lower band) 




7q 



7| 



U 



(B7) 



The degeneracy in the classical (mean-field) dispersion 
along the magnetic zone boundary (et — 0) is lifted by 
the first term in Eq. (B7). The (7r/2, 7r/2) state is pushed 
up more in comparison to the (tt, 0) state, because e^_q 
goes like (singa, -I- sin^j^)^ and {cosq^ — cos qy)'^, respec- 
tively, and therefore the small q contribution dominates 
for k = (7r/2,7r/2). 



2. Ground-state energy 



ent spectral function with identical integrated weight is 
transferred from the upper to the lower band. Focussing 
on the latter, for a state |k©), the first-order correction 
yields, to leading order, 



[5G\'\k,Lj)]BB = [5GY'ik,uj)]AA 



- 1 



(B8) 



and the corresponding incoherent spectral function is 



A incoh 



(u;) 




e2 ) ■ 

k— q/ 



(B9) 

The first-order quantum correction SE^"^ to lower- 
band energy can thus be obtained by simply subtract- 
ing the mean-field contribution corresponding to the (co- 
herent) spectral weight lost from the lower band, and 
adding the energy contribution of the (incoherent) spec- 
tral weight transferred to the lower band. With E^ « 
—A — e^/J7, we obtain 



SE, 



(1) _ 




(BIO) 



Due to mixing of states between the two bands, 
[E^'''' (k, uj)] for a lower-band state yields spectral function 
in the upper band and vice-versa, resulting in transfer of 
spectral weight between the two AF bands. These mod- 
ifications to the lower-band spectral function result in 
a quantum correction to the ground state energy, which 
we evaluate to first order from the corresponding Green's 
function correction. 

While the coherent spectral weight in the lower AF 



band is reduced by ^ J^c 



- 1 



an incohcr- 



Summing over lower-band k states for both spins yields 



for the quantum correction to ground-state energy, which 
is exactly the spin-wave-theory result for the Heisenbcrg 
model. The net lowering in ground-state energy is thus 
a consequence of the incoherent spectral function in the 
lower band, spread over the magnon-energy scale. 
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